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We report results of a search for CPT and Lorentz violation in B°-B° oscillations using inclusive 
dilepton events from 232 million T(4S) — > BB decays recorded by the BABAR detector at the PEP- 
II B Factory at SLAC. We find 2.8<r significance, compatible with no signal, for variations in the 
complex CPT violation parameter z at the Earth's sidereal frequency and extract values for the 



quantities Aa M in the general Lorentz-violating standard-model extension. The spectral powers for 
variations in z over the frequency range 0.26 year - to 2.1 day - are also compatible with no signal. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



It was shown recently [1] that an interacting quantum 
field theory need not be local for CPT violation to imply 
violation of Lorentz invariance. In the general Lorentz- 
violating standard- model extension (SME) [2J, the pa- 
rameter for CPT violation in neutral meson oscillations 
depends on the 4- velocity of the meson [3j . 

We report a search for this effect using T{AS) — * BB 
decays recorded by the BaBAR. detector at the PEP-II 
asymmetric-energy e + e - collider. Any observed CPT 
violation should vary with a period of one sidereal day (~ 
0.99727 solar days) as the T(45) boost direction follows 
the Earth's rotation with respect to the distant stars [J] . 



The physical states of the B"-B° system are 



\b l ) = pVT 
\b h ) = pVT 



z\B°) 

z\B°) 



qVT 

qVT 



z|B°) 
z\B°) 



(1) 



where L (H) labels the "light" ("heavy") eigenstate of 
the effective Hamiltonian. The complex parameter z van- 
ishes if CPT is conserved; T invariance implies \q/p\ = I . 
In the SME, CPT- and Lorentz-violating coupling coef- 
ficients a q j for the two valence quarks in the B° meson are 



contained in quantities Aa. 



n qi 



, a? 2 , where the 



r qi are due to quark-binding and normalization effects. 
The CPT parameter z depends on the meson 4-velocity 
(3^ = 7(1, /3) in each experiment's observer frame as [3[ 



/^Aa, 



»AT/2' 



(2) 



where /3 M Aa M is real and varies with sidereal time due to 
the rotation of j3 relative to the constant vector Aa. The 
magnitude of the decay rate difference AT = Th — ^l 
is known to be small compared to the B°-B° oscillation 
frequency Am = m# — tol ; hence Eq. [2] constrains 



AmRez ~ 2Am(Am/Ar)Imz ~ (3^ Ac 



(3) 



Limits on analogous flavor-dependent Aa M specific to 
K°K° oscillations [j| and to D°D° oscillations @] have 
been reported by the KTeV and FOCUS collaborations, 
respectively. KTeV has also reported a limit on sidereal 

variation of the phase </>_| of the CP- violating amplitude 

ratio ry + _ = A(K L -* n + ir~)/A(K s -> tt+tt - ) 0]. 

We adopt the basis (x, y, z) for the rotating laboratory 
frame and the basis (X , Y, Z) for the Sun-centered non- 
rotating frame containing Aa [8j. Z is parallel to the 
Earth's rotation axis, X(Y) is at right ascension 0°(90°), 
and y is at declination 0° . We take (3^ for each B meson 
to be the T(4S) 4-velocity, and choose z to lie along — (3. 
The event sidereal time t is given by the right ascension 
of z as it precesses around Z at the sidereal frequency 



f2 = 2-7T rad/sidereal-day. We find i = 14.0 sidereal-hours 
at the Unix epoch (00:00:00 UTC, 1 Jan. 1970) from the 
latitude (37.4° N) and longitude (122.2° W) of BaBAR and 
the T(AS) boost ((^7) ~ 0.55 toward 37.8° east of south), 
which also yield cos \ — z ■ Z — 0.628 in Eq. |4| 



/^Acfy = 7 [ Aa - (3 Aa z cos x 

— /3sinx (Aay sin fit + Aax cos Qty\ . 



(4) 



Neutral B mesons from T(4S) decay evolve in orthog- 
onal flavor states until one decays, after which the flavor 
of the other continues to oscillate. We use direct semilep- 
tonic decays (b — » X£v, where t = e or fi) to tag the flavor 
of each B°(B°) by the charge of the lepton £+(£-). The 
decay rate for opposite-sign dilepton (£ + £~) events is 



iv+- 



tx e 



-|A*|/t b0 



{(l + |z| 2 )cosh(AFAt/2) 



(1 



2 ) cos(AmAi) 



(5) 



-2Rezsinh(AFAt/2) + 2Imzsin(AmAi)}. 



We define 1/tb° to be the average neutral B decay rate, 
and At = t + — < - , where t + (t~) is the proper time for 
one of a pair of B mesons to decay to £ + {£r ) . We make 
the approximation sinh(ATAt/2) ~ AFAt/2, which is 
valid for the range |At| < 15 ps used in this analysis. We 
use |AF| = 6 x 10 -3 ps -1 in the cosh(ArAi/2) term, 
consistent with the value reported in Ref. [9j. 

The asymmetry between the decay rates at At > 
and At < compares the probabilities P(B° — ► B°) and 
P(B° — > B°). Omitting second-order terms in z gives 



.4 



CPT 



(At) 



-Re z AT At + 2 Im z sin(AmAt) 
cosh( AFAt/2) + cos(AmAi) 



(6) 



The BABAR detector is described elsewhere [12j . We use 
about 232 million T(45) -> BB decays and 16 fb -1 of off- 
resonance data, from 40 MeV below the T(AS) resonance, 
collected during 1999-2004 to search for variations in z 
with sidereal time of the form 



z = zq + zi cos (fit 



(7) 



For long data-taking periods, any day/night variations in 
detector response tend to cancel over sidereal time. 

We have previously measured [lfj time-integrated val- 
ues of Im z and Re z Ar from the At distribution of the 
same events. Here, we measure Imzo, Rezo Ar, Imzi, 
and Rezi Ar by extending the likelihood fit to include 
the event sidereal time t, and extract values for the 
SME quantities Aa^. In a complementary approach, we 
also measure the spectral power of periodic variations 
in z over a wide frequency band using the periodogram 
method Hj developed to study variable stars. 



The event selection is the same as in Ref. [10| . Briefly, 
we suppress non-BB background by event-shape and 
event-topology requirements, and select events having 
at least two well-identified lepton candidates with mo- 
menta 0.8-2.3 GeV/c in the T(45) rest frame that are 
not part of reconstructed J/ip ,ip(2S) — > e + e~,/i + /i~ de- 
cays or photon conversions. Lepton candidates must have 
at least one ^-coordinate measurement in the silicon ver- 
tex tracker to allow At to be well-measured. We reject 
events in which either of the two highest-momentum lep- 
ton candidates (the dilepton) is classified as a cascade 
lepton from a b — > (c,t) — ► £ transition by a neural- 
network algorithm that uses as input variables the mo- 
menta and opening angle of the two leptons together with 
the event's visible energy and missing momentum. The 
selected dilepton sample comprises 1.18 million opposite- 
sign events and 0.22 million same-sign events. 

We estimate the T(45) decay point in the transverse 
plane with a x 2_n t using the transverse distances to the 
two lepton tracks and the beam-spot. To measure At, 
we assume each lepton originates from a direct B me- 
son decay at the point on the lepton track with the least 
transverse distance to the T(4S). The component Az, 
along the Lorentz boost, of the distance between these 
two points yields At = Az/(j3^)c. For opposite-sign 
events Az — z + — z~; for same-sign events we use |Az|. 

We model the At-distribution of the dilepton sample 
with the probability density functions (PDFs) used in 
Ref. [l0| to represent contributions from B°B° and B + B~ 
decays and non-BB events. The latter are estimated, us- 
ing off-resonance data, to be 3.1% of the sample. The fit 
to data determines that 59% of the BB events are B + B~ 
decays. With minor BB background contributions fixed 
to values from Monte Carlo (MC) simulation, the fit to 
data also determines the fractions of B°B° and B + B~ 
decays that are signal events (~ 80%) with two direct 
leptons, and the fractions (~ 10%) that are events with 
one direct lepton and a&^c^l cascade decay of the 
other B meson. Same-sign dilepton events are retained 
primarily to improve the determination of these fractions. 

Each PDF is a convolution of a decay rate in At with 
a resolution function that is a sum of Gaussians or, for 
events with a cascade lepton, its convolution with one 
or two double-sided exponentials accounting for the life- 
times of intermediate r or Dr s -\ meson states, respectively. 
We use a sum of three Gaussians for signal events. The fit 
to data determines their fractions and also their widths 
except that of the widest, which is fixed to 8ps. For lep- 
tons from different B mesons, our B°B° decay rate con- 
tains z to first-order (cf. Eq. [5]) for opposite-sign events 
and is ex e _|At|/r fl° (cosh(ArAt/2) - cos(AmAt)} for 
same-sign events; for B + B~ decays, it is oc e~' At >' T B ± . 
For leptons from the same B meson, the decay rates are 
exponentials with effective lifetimes determined from MC 
simulation. Dilution factors are included to account for 
wrong flavor tags in cascade decays. 



Each event's timestamp yields the time elapsed since 
the Unix epoch. We use this time, folded over one sidereal 
day and shifted in phase by 14.0 sidereal-hours, for t. 

We extract z from a two-dimensional maximum likeli- 
hood fit to the opposite-sign and same-sign data events 
binned separately in At and t. The likelihood function in 
At for each of the 24 sidereal-time slices contains a com- 
mon sum of the PDFs, and z varies with t as in Eq. [7] 
The likelihood fit corresponds to Aqpt m Eq. [6l We 
obtain the values for z and 4> reported in Table U (up- 
per left). The statistical correlation between Imzo and 
Rez Ar is 76%; between Imzi and Rezi Ar it is 79%. 

Table|T]shows the sources of systematic uncertainties in 
the asymmetry parameters. Separate contributions are 
added in quadrature in the totals. We vary separately 
t b o, Tg±, and Am by la from their known values [13j , 
and vary |Ar| over the range 0-0. lps -1 to allow 3a 
deviations from the value reported in Ref. [9| . Fixed pa- 
rameters in the PDF resolution functions for non-signal 
events are varied separately by 10%, motivated by a com- 
parison of resolution parameters fitted to signal events in 
data and MC simulation. The fractions of the Dt s \ meson 
components in background cascade decays are also varied 
by 10%. The effects of possible internal misalignments 
of the silicon vertex tracker (SVT) and uncertainty in 
the absolute z-scale are evaluated in B°B° MC samples. 
The clock that sets the event timestamps is governed by 
the PEP-II master oscillator, which is stable to within 
0.001% of its set frequency. Resynchronization of the 
clock with U.S. time standards at intervals of less than 
four months limits relative sidereal phase errors to less 
than 0.2%. Another small uncertainty in sidereal phase 
arises in calculating the Y(4S) boost's right ascension. 
We use e + e _ — ► n + fi~('j) data events, with true Az = 0, 
to check for sidereal variations in measured Az that could 
mimic a Lorentz- violation signal. The measured ampli- 
tude (0.022 ± 0.025) ^m and mean (0.030 ± 0.018) (im are 
sources of negligible uncertainties. At the solar-day fre- 
quency, the amplitude is (0.028 ± 0.025) /xm. 

In Fig. Q] we plot the sidereal-time dependence of the 
measured asymmetry A™ c j^f, for the opposite-sign dilep- 
ton events with |At| > 3ps, thereby omitting highly- 
populated bins where any asymmetry is predicted to be 
small. Figure [2] shows confidence level contours for Imzi 
and Rezi Ar. The significance for sidereal variations in 
z, characteristic of CPT and Lorentz violation, is 2.8er. 

The results of the fit described above are compatible 
with the SME constraint Rez Ar ~ 2AmImz (Eq.^J) for 
Am = 0.507 ps _1 [Hj]. We repeat the likelihood fit sub- 
ject to this constraint. The asymmetry in Eq.[6]becomes 



.4. 



CPT 



(At) 



2 Imz {-AmAt + sin(AmAi)} 
cosh(ArAt/2) + cos(AmAt) 



(8) 



We obtain the results reported in Table I (right). The 
statistical correlation between Imzi and <j> is 48%. The 
significance for sidereal variations in z is again 2.8<r. We 



TABLE I: Asymmetry parameter values, with statistical errors, for Acpt in Eq. [6] (upper left) and with SME constraint in 
Eq. [8] (upper right). Equation [7] implies zi — > — zi for — > <^> + 7r. Systematic uncertainties are shown in lower part of Table. 







Without SME constraint 




With SME constraint 


Acpt parameter- 


Imz 


Re zo Ar 


Imzi 


Rezi Ar 





Imz 


Imzi 







(xl(T 3 ) 


(xl0 -3 ps -1 ) 


(xlO -3 ) 


(xl0 -3 ps -1 ) 


(rad) 


(xlO -3 ) 
-5.2 ±3.6 


(xlO -3 ) 
-17.0 ±5.8 


(rad) 


Value from fit 


-14.2 ±7.3 


-7.3 ±4.1 


-24 ±11 


-18.5 ±5.6 


2.63 ±0.31 


2.56 ±0.36 


Systematic effects 


















t b° i t b± > Am, af 


±0.7 


±0.4 


±0.6 


±0.5 


±0.05 


±0.4 


±0.7 


±0.01 


SVT alignment, z scale 


±0.6 


±1.5 


±2.0 


±1.1 


±0.20 


±1.7 


±1.4 


±0.15 


PDF resolution models 


±2.0 


±1.0 


±2.5 


±1.2 


±0.02 


±0.8 


±1.0 


±0.01 


Background fractions 


±0.1 


±0.1 


±0.2 


±0.2 


±0.01 


±0.2 


±0.3 


±0.01 


Sidereal phase 


±0.0 


±0.0 


±0.0 


±0.0 


±0.03 


±0.0 


±0.0 


±0.03 


Total syst. error 


±2.2 


±1.8 


±3.3 


±1.7 


±0.21 


±1.9 


±1.9 


±0.15 



< 0.03 ; 
0.02; 
0.01 

-0.01 
-0.02 ; 

-0.03 



I. .. i ... i ... i ... i ... i ... i ... i ... i ... i ... i ... i . . 
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FIG. 1: Asymmetry Aq c p T for opposite-sign dilepton events 
with I At I > 3ps versus sidereal time. The sample includes 
event types, e.g. B + B~ decays, for which Acpt — 0. The 
curve is a projection, for \At\ > 3ps, using results of the 
two-dimensional likelihood fit for |Ai| < 15 ps. 
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FIG. 2: Contours indicating la, 2a, and 3a significance, 
around the central values of Imzi and Rezi Ar (solid circle). 



obtain consistent results for Imzo, Imzi, and <f> when 
second-order terms (Eq.fSJ) of form |z| 2 = p 2 cos 2 (Sit + (f>) , 
motivated by finding |Imzi| > |Imz |, are included in the 
likelihood fit to data with p 2 as a free parameter. 

We use Eqs. |3l HI and [7J to extract the SME quantities 



Aa 



0.30Aaz 
Aa x 
Aay 



(-3.0 ±2.4) (Am/ Ar) x KT 15 GeV, 
(-22±7)(Am/Ar) x KT 15 GeV, 
(-14±i§)(Am/AT) x 10 -15 GeV. 



We now use the periodogram method [ll| to compare 
the spectral power for variations in z at the sidereal fre- 
quency with those in a wide band of surrounding frequen- 
cies. The spectral power at a test frequency v is 



P(v) 



1 



Nal 



N 



"li-KvTi 



(9) 



where the data, comprising N measurements Wj made at 
times Tj, have variance a\,. Here, Tj is the time elapsed 
since the Unix epoch for opposite-sign dilepton event j, 
and the weights Wj — AmAtj — sin(A?7iAtj) are suited to 
the study of periodic variations in z according to Eq. [5J 
In the absence of an oscillatory signal, the probabil- 
ity that P(v) exceeds a value S at a given frequency is 
exp(— S); if M independent frequencies are tested, the 
largest P(v) value exceeds S with probability 



Pr{P max (v)>S;M} = l-[l- 



-,-s 



M 



(10) 



We use 20994 test frequencies from 0.26 year 



to 

2.1 solar-day -1 , spaced by 10~ 4 solar-day -1 . This over- 
samples the frequency range by a factor of about 2.2 and 
avoids underestimating the spectral power of a signal. 
The number of independent frequencies is about 9500. 

Figure[3]shows the periodogram we obtain. The largest 
spectral power is i- > ma x(i / ) = 8.78, for the test frequency 
v = 0.46312 solar-day -1 . With no signal, the prob- 
ability of finding a larger spectral power in our peri- 
odogram is 76%. Interpolation to the sidereal frequency 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
Frequency (1 /solar-day) 

FIG. 3: Periodogram for opposite-sign dilepton events. The 
solar-day and sidereal-day frequencies are indicated by the 
left and right triangles, respectively, in the inset. 



(~ 1.00274 solar-day^ 1 ) yields P[y) = 5.28, a value that 
is exceeded at 78 test frequencies. At the solar-day fre- 
quency, where any effects due to day/night variations in 
detector response should appear, P(v) = 1.47. 

In conclusion, we report results of a search for side- 
real variations in the CPT violation parameter z that 
are complementary to our previous time- integrated mea- 
surements [lCf using the same events. Neither the like- 
lihood fits nor the periodogram method detect asym- 
metries large enough to provide evidence for CPT and 
Lorentz violation. We have constrained the quantities 
Aa M of the Lorentz-violating standard-model extension 
that parameterize CPT violation in B a -B° oscillations. 
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